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ABSTRACT 



Context. High-energy y-rays propagating in the intergalactic medium can interact with background infrared photons to produce e + e~~ pairs, 
resulting in the absorption of the intrinsic y-ray spectrum. TeV observations of the distant blazar 1ES 1 101-232 were thus recently used to put 
an upper limit on the infrared extragalactic background light density. 

Aims. The created pairs can upscatter background photons to high energies, which in turn may pair produce, thereby initiating a cascade. The 
pairs diffuse on the extragalactic magnetic field (EMF) and cascade emission has been suggested as a means for measuring its intensity. Limits 
on the IR background and EMF are reconsidered taking into account cascade emissions. 

Methods. The cascade equations are solved numerically. Assuming a power-law intrinsic spectrum, the observed 100 MeV - 100 TeV spectrum 
is found as a function of the intrinsic spectral index and the intensity of the EMF. 

Results. Cascades emit mainly at or below 100 GeV. The observed TeV spectrum appears softer than for pure absorption when cascade 
emission is taken into account. The upper limit on the IR photon background is found to be robust. Inversely, the intrinsic spectra needed to 
fit the TeV data are uncomfortably hard when cascade emission makes a significant contribution to the observed spectrum. An EMF intensity 
around 10~ 8 nG leads to a characteristic spectral hump in the GLAST band. Higher EMF intensities divert the pairs away from the line-of-sight 
and the cascade contribution to the spectrum becomes negligible. 
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1. Introduction 

The observed very high energy (VHE) spectra of extragalac- 
tic sources are attenuated by pair production (PP) on back- 
ground photon fields. At energies e = 100 TeV, interactions 
with CMB photons make the universe opaque beyond a few 
Mpc dGould & Schrederlll966t) . In the 1-10 TeV range, the tar- 
get photon field is the infrared extragalactic background light 
(EBL), with an undetermined horizon z ^ 0.2 due to un- 
certainties in the EBL density at optical to UV wavelengths. 
Conversely, observations of absorbed VHE spectra of blazars 
can constrain the EBL density at t hese wavelengths, provided 



their intrinsic emission is known (Stecker et al. 1992; Biller 



1995). Recently, the HESS collaboration used its observations 
of 1ES 1101-232 (z = 0.186), together with a reasonable as- 
sumption on the intrin sic spectrum, to estim ate a stringent up- 
per limit to the EBL dAharonian et al. I l2006l) . This estimation 
did not consider emission from the cascade initiated when the 
created pairs upscatter EBL photons back t o VHE energies 



dProtheroel 1 19861: IProtheroe & Stanevl 1 19931: Aharonian et al 
Il994t lBillerlll995t lAharonian et al.ll2002h . Cascade emission 
may make the universe appear more transparent than under 



the assumption of pure absorption. Inversely, including cascade 
emission when deconvolving for propagation effects on a given 
EBL leads to intrinsic spectra different from the pure absorp- 
tion case. 

The impact of this emission on the EBL upper limit, 
as derived from the 1ES 1101-232 spectrum, is considered 
here. Electrons may diffuse on the extragalactic magnetic 
field ( EMF), causing their emissions to be lo s t for the ob - 
server (Protheroe 1986; A haronian et al.1 1 19941 IPlagal Il995h . 
The EBL and EMF governing the propagation of the cas- 
cade are described in §2. The cascade equations and numer- 
ical method are described in $3] and applied to the case of 
1ES 1101-232 in §4. The implications on the EBL and EMF 
limits are set out in §5. 

2. Extragalactic backgrounds 

In addition to the CMB, the diffuse photon background is 
constituted of integrated emission from stars (peaking around 
2//m, see Fig. [TJ and heated dust (peaking around 200 fim). 
The EBL spectral energy distribution is difficult to measure di- 
rectly because of the strong zodiacal and galactic foregrounds. 
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Fig. 1. EBL+CMB photon density (at z-Q) used in this work. 
The HESS 2006 and Prima ck 2005 EBL derive from a simu- 
lation of galaxy formation (IPrimack et al. 1999, I2005L respec- 
tively ). The Spitzer 200 6 EBL is a best fit to available observa- 
tions (Dol e et al. 2006, from which the measurements shown 
here were also taken). 
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Fig. 2. Mean free path A PP for PP on the EBL+CMB back- 
ground as a function of the VHE photon energy. The mean 
free path A\c and IC cooling length Cic of the pairs on the 
CMB is also shown as function of electron energy. Cic flat- 
tens (then rises) at high energies as the interaction enters the 
Klein-Nishina regime. The electron then loses its energy in a 
single interaction (Aic ~ Cic). IC losses of the pairs take place 
on a small scale compared to the y-ray attenuation length for 
photon energies < 300 TeV. The electron gyroradius R L for var- 
ious EMF intensities is indicated by dashed diagonal lines. The 
pairs are expected to be isotropised by the EMF for energies 
and B intensities such that < Cic- 



Lower limits have been set using source counts while SED 
shap es have been derived from simulations of galaxy fo rma- 
tion JPrimackll200l lLagache et aTll2003l: IXu et alJl200lh. The 
EBL s hape HESS 2006 (Fig.UJ was computed bvlPrimack et alJ 



1999) and normalised by a factor 0.45 in lAharonian et al.l 
20061) to account for the TeV observations of 1ES 1101-232. 
Primack 2005 refers to an EBL normalised instead to fit the 
lower limit set by galaxy counts. The Spitzer observations sug - 
gest higher fluxes in the 10 - 1000 yum range ( Dole et al.l20 06). 
This affects attenuation above 20 TeV but has been verified to 
have no consequence on the results presented here. 



The created pairs can be deflected from the line-of-sight 
by an extragalactic magnetic field (EMF). Faraday rotation 
and synchrotron emission in radio yield estimates of magnetic 
fields in galaxies (roughly > 10 nG), or in c lusters (< 0.1 - 
1 nG) and even some super - clusters (< nG) dKronbergl 1994; 
WidrowlEool I Valleebooi) . The EMF outside these structures 
is un constrained and may be as low as 10~ 19 nG (iFan et al 



2003, and references therein). For such very weak EMFs, the 
deflection of electrons due to IC interactions is negligible and 
the cascade occurs along the l ine-of-sight with a short de- 
lay of the seco ndary emission dPlagal 119951: ICheng & Chengl 
1996k iDai etalT] 2002). Diffusion on a stronger EMF creates a 
halo around y-ray sourc es and isotropizes the cascade emission 



(lAharonian et alj|1994l) . This occurs when the gyroradius 
of the pairs is much lower than their Compton cooling length 
Cic = E(dE/dl)7}, Since mostly CMB photons are upscat- 
tered, the minimum B required to isotropise pairs of energy E is 
3 lO -6 ^ v (l +z) A nG. Much of the isotropic re-emission is lost 
to the observer and the pairs diffuse on a scale ~ (R^Cic) 1 ^ 2 - 
For intermediate EMFs, the TeV electrons in the beamed rela- 
tivistic jet are deflected by ~ C\c/Rl- Halo sizes ^0.1° could 
be resolv ed by y-ray detectors and u sed to estimate the EMF 
intensity ( Neronov & Semikozl 20061) . Photons in 0.1° haloes 
have propagation ti mes varying by ~ 10 5 years, averaging out 
any time variability dFan et alj |2003). In the following, the cas- 



cade emission is assumed to be unresolved from the source and 
delays are not considered. The TeV emission detected by HESS 
from 1ES 1101-232 appears to be at a low flux level with no 
significant variability. 

3. Cascade equations 

The cascade is described by a set of two coupled equations 
involving the photon energy density rip(s) and the electron 
(positron) energy density rtE(E): 



cd,np 



i r 

-—rip + CB 

App Je 

1 



cd t n E — ri£ + 2 

^ic Je 



Gic(e,e)n E (e)de 



Gpp(e, E) n e (e)de 



r 



(i) 



(2) 



Gic(e, e — E) n E {e)Ae 



The first term in both equations is the sink term due to 
PP (Eq. [TJ or IC losses (Eq. |2j. App and A\c are the mean 
free path for each interaction. The second term is the source 
term corresponding to cascade emission (Eq. [TJ or pair cre- 
ation (Eq. [2] with a factor 2 for the pair). The cascade 
emission factor cb is 1 when the EMF is ignored, and ap- 
proximated to when the electron population is considered 
isotropised. The pair production term is written in terms of 
Gpp(e, E) = J dE<rpp(s, £b)m(sb)^£b, where Sectpp is the dif- 
ferential cross-section and u is the photon background energy 
density (EBL+CMB). The IC radiation term Gic(£', s) is de- 
fined similarly. The third term in Eq. [2] reflects IC cooling of 
electrons from higher energies. All of these terms are functions 
of z- 

The integrated cross-s ections for PP and IC on iso tro pic tar- 
get photons are taken from Gould & Schredei ( 19661) and Jones 
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Fig. 3. 1ES 1101-232 observed and modelled spectra with 
a maximal (HESS 20006) EBL and including cascades with 
no magnetic field. HE SS observations points are in black 
dAharonian et al]|2006l) . Markers indicate the attenuation only 
observed spectrum (circles) and the corresponding intrinsic 
spectrum (crosses), whereas the lines indicate the observed 
spectra with cascade emissions. Intrinsic spectra are in the form 
of vF v cc E° dv and adjusted to the data. Cascade emission 
accumulates at 100 GeV and below, softening the spectra com- 
pared to pure absorption. The HESS upper limit on the EBL 
remains valid after taking the full emission from cascades into 
account. 



( 1967|). An alytic expr e ssions of the differential cross- sections 
derived by IZdziarskil dl988l) for background densities in the 
form of blackbodies or power laws are used to calculate Gpp 
and Gic- The cascade equations are solved numerically by com- 
bining np and he into a single vector V defined on a logarithmic 
scale of energies (eo<T'X from en=10 7 eV up to 10 17 eV in 250 
steps (Thus £ = (10 17 /10 7 ) I/2S0 ). To ensure energy conserva- 
tion, the integrals on Gpp and Gic are calculated as 



Gic(e, s) n E (e)de 



V Jr 1 ' 2 



eo£ k uG\c(eo£ k u, s)du 



£1/2 



-1/2 



(3) 

The cascade equations may then be rewritten as a matrix P 
acting on the vector V : V(t + 6t) = exp(5fP)V(f) (exp is de- 
veloped to the 4 th order in 5t). The terms in P are of the or- 
der of or less, hence it is enough to take steps of size 
c6t = 0.1 kpc, updating the matrix P(z) every Sz = 0.001 with 
d z = H n (l +z)[Q M (l+z) 3 + QA+(l-QM-QA)( l+z) 2 l 1/2 df 
and v alues for Hq,Q.m and £2 a taken from WMAP (iPeiris et al 
2003b . Thus, at z = 0.2, Sz ~ 3 10 4 aJf. 



4. Application to 1ES 1101-232 

The SED of the attenuating EBL can be deconvolved from y- 
ray observations of extragalactic sources (TeV blazars), given 
a priori knowledge on the intrinsic spectra. Modelling ob- 
served spectra as power-laws, the effect of PP is to soften the 
intrinsic spectral index, increasingly so with EBL intensity. 
Hence, using observations of the farthest TeV blazar and as- 
suming the hardest possible intrinsic spectrum puts an upper 
limit on the EBL responsible for attenuation. Current theoreti- 
cal understanding of shock acceleration limits the intrinsic par- 
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Fig. 4. 1ES 1101-232 observed and modelled spectra with a 
minimal (Primack 2005) EBL and including cascades with no 
magnetic field. The intrinsic spectrum is now adjusted to the 
HESS data leaving the spectral index F free. For pure absorp- 
tion, the best index is F — 1.95. With cascades, the index hard- 
ens as IC emission softens the propagated spectrum. For high 
€m, the best index softens again so as to limit the amount of 
cascading but the fit worsens. Significant cascading on the min- 
imal EBL and in a very weak EMF implies a very hard input 
spectrum in order to account for the observations. 



tide distribution in blazars to a power-law of index no harder 
than a 1.5 and correspond ingly, an intrinsic pho ton spectrum 
dN oc E r dE with Y > 1.5 dAharonian et alj|2006l) . 

1ES 1101-232, atz = 0.186, is currently the farthest known 
TeV source and was used by the HESS collaboration to set an 
upper limit to the EBL corresponding to the HESS 2006 SED 
shown in Fig. 1. The comparison between a T=1.5 power-law 
attenuated by the HESS 2006 EBL (without cascade, c^=0) 
and the data is shown in Fig. 3, reproducing the results of 
Aharonian et al. (2006). Attenuated spectra taking into account 
the full cascade emission with c B -\ (i.e. a null EMF) are also 
shown for various values of the maximum energy e M to which 
the intrinsic power-law extends. Since cascades initiated at 
higher energies increase the photon populations in lower ones, 
one might expect the final spectra to appear harder than for pure 
absorption. However, because IC occurs predominantly on the 
CMB, the cascade emission accumulates below 100 GeV, soft- 
ening the spectrum between 100 GeV and 1 TeV. High values 
of £m lead to more cascading and more softening. The^- 2 values 
suggest 6m < 15 TeV, although further observations, particu- 
larly above 1 TeV, would be necessary in order to confirm this. 
For such low €m values, not many photons initiate cascades. 
For higher e^, the softening is such that a lower EBL would 
be neede d to match the data. Th us the HESS 2006 upper limit 
found by I Aharonian et al. (2006) holds strong, even in this ex- 
treme limit where all the cascade emission is received by the 
observer. 

Inversely, the intrinsic y-ray spectrum at the source can 
be obtained given some assumption on the intervening EBL. 
Using the lower limit on the EBL set by galaxy counts (Primack 
2005 in Fig. 1) gives a limit on how soft the intrinsic spectrum 
can be. For pure absorption, the best fit has Y - 1.95 + 0.19 
(Fig. 4). As expected, this is softer than the F = 1.5 assumed 
above, yet still suggests that a good fraction of the y-ray energy 
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Fig. 5. Observed spectra for the Primack 2005 EBL and various 
EMFs between 1CT 9 and 10~ 6 nG. The spectra are adjusted to 
the HESS points, leaving F free but fixing £m = 10 TeV (top) 
or 20 TeV (bottom). For high EMFs the cascade emission does 
not reach the observer and the spectrum is identical to the pure 
absorption case. For very low EMFs the results are similar to 
those shown in Fig. 4. Intermediate values lead to a more or 
less pronounced bump at 1-100 GeV energies over the intrin- 
sic continuum. Hatches represent projected 5-cr sensibilities for 
one year of observation with GLAST (< 30 GeV) and 50 hours 
with HESS-2 (> 50 GeV). 



in 1ES 1 101-232 is output above a TeV. A hard T < 2 intrinsic 
spectrum is needed if cascade emissio n is to contribute signif - 
icantly to the low-energy continuum dAharonian et alj|2002l) . 
1ES 1101-232 is the first blazar where the intrinsic spectrum 
is constrained to be hard enough for this, even in the minimal 
EBL limit. 

Including cascade emission in the fit (Fig. 4) hardens even 
more the intrinsic spectrum as the cutoff e M increases and cas- 
cades contribute more and more to the observed spectrum. For 
higher eu, the best fit T increases again to mitigate the pro- 
nounced softening from the strong cascading but the fit wors- 
ens. This also holds for (implausibly) high values of €m > 
100 TeV, for which cascade emission largely dominates at a 
few TeV. The hard intrinsic spectra found here, assuming the 
Primack 2005 is indeed the minimum possible EBL, suggest 
either that 6m is not greater than a few TeV, so that there is little 
cascade emission in the TeV range, or that a large part of the 
cascade emission is lost due to diffusion on the EMF 

As discussed in §2, the electron diffusion on the EMF de- 
pends on the ratio R^/Cic- The effect on the observed spectra 
is now taken into account by setting c B =0 when 7? L /Cic<300 



(corresponding to a maximum deviation on the line-of-sight 
of 0.1°-0.2° equal to the best GLAST angular resolution) and 
cb-1 otherwise. For example, an EMF of 10~ 6 nG means that 
emission from electrons of energy E $ 20 TeV is suppressed. 
This will lead to low-energy cutoff in the cascade spectrum as 
only emission from pairs above a certain energy reaches the 
observer. The overall spectrum appears as a hump between 
y 2 hvcMB (with y the Lorentz factor of the electrons for which 
Rl = 300Cic) and 100 GeV (above which absorption dom- 
inates). Hence, a non-zero EMF leads to a reduction of the 
overall cascade emission seen by the observer (compared to 
Figs. 3-4) but can also lead to a well-defined signature above 
the continuum. 

Figure shows the observed spectra for a Primack 2005 
EBL and for EMF intensities between 10~ 9 and 10~ 6 nG. The 
intrinsic power-law index was left free but its cutoff £m was 
fixed at either 10 TeV or 20 TeV. The best fit index F is then 
found for each value of the EMF. In both cases, the spectra 
for an EMF <; 10~ 6 nG are not much different from the pure 
absorption case as most of the cascade emission is isotropised 
and lost to the observer. With 6m=10 TeV, the best-fit intrinsic 
slopes are flat in vF v and the cascade emission is essentially 
indistinguishable from the GeV continuum for any value of the 
EMF. The intrinsic emission is assumed here to be a simple 
power-law over the whole energy range. More realistic mod- 
elling would result in a curved intrinsic Compton component. 
The cascade emission might then be more readily identifiable 
over an intrinsic continuum rising from GeV to TeV energies. 

Stronger cascading, as a result of a higher cutoff energy 
€m and/or a higher EBL density, makes the hump apparent 
for the same reason. The intrinsic spectrum is then necessar- 
ily much harder, enabling the contribution from the cascade to 
stand out over the continuum. The bottom panel of Fig. 5 shows 
that EMF intensities of 10 9 — 10 8 nG can be identified using 
GLAST and HESS-2 if 6m=20 TeV. Cascade emission is not 
diluted for EMF intensities weaker than 10~ 9 nG and there is 
no spectral feature to measure the EMF. Surprisingly, in most 
cases 1ES 1101-232 is only slightly above the GLAST one- 
year detection limit. Unless they become active and flaring, low 
flux state blazars detected by HESS such as 1ES 1 101-232 are 
likely to be difficult to detect with GLAST, illustrating the ad- 
vantage provided by the large collecting area of ground-based 
Cherenkov arrays (but at higher energy thresholds). Similar 
results are obtained by keeping €m at 10 TeV but using the 
stronger HESS 2006 EBL. However, in this case, the fitted in- 
trinsic slopes are very hard (F * 1.1) when the EMF intensities 
are lower than 10~ 7 nG. 

The softest values of F, which are the most plausible given 
the present knowledge on blazars, favour values of the EMF 
higher than 10~ 6 nG and/or a cutoff energy below 20 TeV. VHE 
emission from nearby, little-attenuated blazars can be investi- 
gated for evidence of cutoffs at energies > 20 TeV — although 
it should be noted that e.g. HESS observations of Mkn 421 
(z = 0.03) taken at a high flux actually me asure an exponen- 
tial cutoff at 3 TeV (lAharonian et al.l 120051) . EMF intensities 
^ 10~ 6 nG are consistent with measures inside clusters and 
super-clusters. Such structures may reach 10-50 Mpc in size, 
which is greater than the attenuation length for y-rays above 
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50 TeV. Furthermore, the largest voids, where the EMF is ex 
pecte d to be very small, have a size (20 h~ x Mpc. iPatiri et al 



20061) . smaller than the distance to 1ES 1 101-232. Hence, cas- 
cades are likely to be initiated inside walls. As Cic is only of 
the order of 1 Mpc, such cascades reemit most of their energy 
within the confines of the clusters, and thus are subject to dif- 
fusion. In this case, the cascade emission can only be detected 
by resolving the faint halo surrounding the y-ray source. 

5. Conclusion 

The impact of extragalactic cascade emissions on the GeV-TeV 
spectrum of 1ES 1101-232 has been investigated and shown 
to soften the observed spectrum in the TeV range compared to 
pure absorption. This occurs because most of the cascade emis- 
sions occurs at 100 GeV and below. As a result, the upper limits 
on the EBL determined by HESS are strengthened in the sense 
that taking cascades into account would lead to harder intrinsic 
spectra than judged plausible, or to a reduced EBL upper limit. 
Inversely, using lower limits on the EBL coming from galaxy 
counts, the intrinsic spectrum of 1ES 1 101-232 is found to have 
T < 1.95, with very hard values if there is an important con- 
tribution from cascade emission. This is at odds with current 
theoretical and observational understanding of blazars. A cut- 
off S= 10 TeV in the intrinsic spectrum would limit the cascade 
contribution. This contribution would also be quenched if the 
EMF intensity is greater than 10~ 6 nG, as expected away from 
voids. A lower EMF increases the amount of cascade emission 
reaching the observer in the GeV band, with a signature in the 
GLAST band for intensities ~ 10~ 8 nG — but at the price of a 
hard intrinsic spectrum so as to fit the HESS observations. 
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